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ABSTRACT 

We present a new analysis of the physical conditions in three Herbig-Haro complexes (HH 54, 
HH 212, and the L 1157 protostellar jet) using archival data from the Infrared Array Camera 
(IRAC) on the Spitzer Space Telescope. As described in detail in Paper I, the emission observed 
using the 4.5-/Ltm filter is enhanced in molecular shocks (T=1000-4000 K) at relatively high 
temperature or densities compared with that observed with the 8.0-^m filter. Using these data 
sets, we investigate different distributions of gas between high and low temperatures/densities. 
Our analysis reveals the presence of a number of warm/dense knots, most of which appear to 
be associated with working surfaces such as the head of bow shocks and cometary features, and 
reverse shocks in the ejecta. These are distributed not only along the jet axis, as expected, but 
also across it. While some knotty or fragmenting structures can be explained by instabilities in 
shocked flows, others can be more simply explained by the scenario that the mass ejection source 
acts as a "shot gun", periodically ejecting bullets of material along similar but not identical 
trajectories. Such an explanation challenges to some degree the present paradigm for jet flows 
associated with low-mass protostars. It also give clues to reconciling our understanding of the 
mass ejection mechanism in high and low mass protostars and evolved stars. 



Subject headings: ISM: Herbig-Haro objects 

1. Introduction 

Collimated jets have been observed in a num- 
ber of low-mass protostars over a variety of stellar 
masses and evolutionary stages. These jets have 
traditionally been seen via the thermal excitation 
of atoms and ions in partially ionized shocked re- 
gions, i.e., line emission at optical wavelengths 



ISM: jets and outflows — infrared: ISM 

The popular magneto-centrifugal models beau- 
tifully explain the driving and collimation of these 



jets (Shu et al. 2000 Konigl and Pudritz 20001 



(Ha, [Oi] etc., see Bally et aL||2007 for a review). 
More recently, jets associated with the youngest 
protostars have been observed in molecular emis- 
sion at radio (e.g., CO, SiO, SO) and infrared 
wavelengths (e.g., H 2 ) (see Arce et aL"||2007 Bally 
et al.||2007l for reviews). 
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Magnetic fields are assumed to be coupled either 
to the star and star-disk interface (in the so-called 
X-wind model) or the disk surface (in the disk 
wind theory). At the inner edge or surface of 
the accretion disk, the gas is lifted from the disk 
and accelerated due to a combination of centrifu- 
gal force and magnetic fields aligned with the disk. 
Such a flow will spin, and at larger distances the 
spinning motion will twist the magnetic field, pro- 
viding a "hoop stress" which collimates the out- 
flowing gas. Possible non-steady alternative mech- 
anisms for the driving include magnetic pressure 
(e-g- 



Uchida and Shibata 


1985 


Machida et al. 


and stress (e.g., 


Hayashi et al.||1996 Good- 



son et al.|1999[ ) . Although the driving mechanism 



is not yet clear, astronomers are reaching the con- 
sensus that jets are magneto- hydrodynamically 
driven and powered by mass accretion. 
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In their simplest form some of the above the- 
ories predict that the ejecta of the jet are dis- 
tributed as a continuous flow or, as seems more 
likely due to unsteady mass accretion, a chain 
of knotty structures moving in opposite direc- 
tions. In reality, the observations of protostcllar 
jets show a much more complex variety of features 
(knots, bow shocks, bubbles etc.) than predicted 
by idealized models. These features are usually at- 
tributed to hydrodynamic/magnetohydrodynamic 
instabilities of the flows, or time variability of the 



ejection in mass, velocity and/or direction (Bally 
et al.|2007||Arce et al.|2007[ ). Mass ejection that is 



variable in mass and velocity can be attributed to 



instabilities in the accretion disk (e.g., Bell and Lin 
1994 Balbus and Hawley|19 98) or the episodic re- 



connection of magnetic field lines between the pro- 
tostar and disk (e.g., Hayashi et al.|1996 Goodson 
etlL|[l999| |. 



While low-mass protostars tend to exhibit colli- 
mated jets, high-mass protostars tend to show out- 
flows with wide opening angles ( Arce et aT||2007 |. 
Studies of warm shocked gas show that such out- 
flows, particularly the most energetic ones, are 
highly turbulent. The best known example is the 
OMC-1 outflow, which consists of many chains of 
knots (finger structures) which have been ejected 



into the ambient cloud (Allen and Burton 1993 



Nissen et al. 2007). In contrast, the Cepheus A 
outflow complex is associated with chains of bow 
shocks in near-infrared H2 emission, and these can 
be explained with a pulsed, precessing jet (Cun- 
ningham et al.|2009 ). It is likely that outflows as- 



sociated with high-mass protostars are also pow- 
ered by mass accretion, however, it is not clear 
how their driving mechanism is related to that of 
the collimated jets seen toward low-mass proto- 
stars ( |Arce et al.||2007[ ). 

In general, observational studies of the driving 
of protostellar jets and outflows have been ham- 
pered by several issues. Their driving engines are 
expected to be too compact to be directly resolved 
without interferometers (i.e., they are within 2-3 
AU of the protostar, corresponding to 15-20 mil- 
liarcsec in the nearest star forming regions at 140 
pc — |Shu et al.||2000l |Konigl and Pudritz|[2000 



Machida et al. 2008; Hayas hi et al.|1 996 Goodson 



et al.||1999 ). Moreover the "engine" is often heav- 
ily embedded, limiting the observing wavelengths 
to radio and far infrared. Finally, in many cases 



we are not directly observing the ejected gas, but 
rather the gas that is shocked or compressed due 
to interactions between the jet and ambient gas 
( pally et al.|2007| |Arce eTaT1[2007l ). 



The InfraRed Array Camera (IRAC) on the 
Spitzer Space Telescope has been used to study 
a variety of protostellar jets and outflows (e.g., 



Noriega-Crespo et al. 2004| |Smith et al. |2006 


Walawender et al.||2006 


|Davis et al.||2007| |Qiu 


et al.||2008| Ybarra and 


Lada||2009| Neufeld et al. 


2009). Takami et al. ( 


2010 hereafter Paper I) 



investigated the nature of infrared emission with 
IRAC in six Herbig-Haro complexes. In Paper I 
we showed that thermal H 2 emission at T=1000- 
4000 K can explain the morphological similarities 
and differences observed between the four bands, 
and also the spectral-energy distributions (SEDs) 
observed in some regions. SEDs in other regions 
show excess emission at 4.5 /an, presumably be- 
cause of excitation of the CO band at high densi- 



ties. See also 


Smith and Rosen (2005|); 


Ybarra and 


Lada 


(2009); 


Neufeld and Yuan (2008 


1 for similar 



studies. The observations and LTE/non-LTE cal- 
culations of thermal H 2 emission indicate that the 
emission at 3.6 and 4.5 /im tends to be associ- 
ated with regions at high temperatures or densi- 
ties, while that at 8.0 /zm tends to be associated 
with regions at low temperatures or densities. Pa- 
per I also suggests that the /4.5^m/^3.6/im ratio 
is useful for discriminating between Herbig-Haro 
knots and foreground/background stars: this ra- 
tio is, >1.6 for Herbig-Haro knots; ~0.6 for fore- 
ground stars with Ak=0; ~1.1 for background 
stars with Ak=5. 

In this paper we extend our analysis for three 
objects (the L 1157 protostellar jet, HH 54, and 
HH 212) to understand their interaction with the 
ambient gas and investigate the nature of their 
mass ejection. We analyze the distribution of 
shocked gas at high temperatures/densities using 
the differing spatial distributions of emission at 
4.5 and 8.0 /mi. The remaining part of the pa- 
per is organized as follows. In §2 we summarize 
the data set we have used. In §3 we describe our 
approach for discriminating between warm/dense 
and cool/diffuse regions, and show how these are 
distributed in the individual objects. In §4 we dis- 
cuss their physical conditions based on the IRAC 
colors, and also the spectra obtained using the 
Spitzer InfRared Spectrograph (IRS) for the L 
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1157 jet (Neufeld et al. 2009 Nisini et al. 20101 



In §5 we discuss the origin of these structures, and 
implications for the distribution of the ejecta from 
the source and their driving mechanism. In §6 we 
present our conclusions. 

2. Data 

Archival data in four IRAC bands (3.6, 4.5, 
5.8, and 8.0 /xm) were obtained for the L 1157 jet 
(d -25 pc, |Looney et al.||2007]) , HH 5 4 (d -200 



Caratti o Garatti et al.| 



Lee et al. 2008 



2006), and HH 
and refer- 



pc, see 

212 (d -400 pc 
ences therein). All the data had been reduced 
with the post-BCD pipelines developed by IPAC. 
The mean FWHMs of the point response functions 
(PRFs) are 1.66, 1.72, 1.88, and 1.98" for the four 
bands, respectively. Previous publications with 



these data include Looney et al. (2007); Neufeld 



et al.| (2009) for the L 1157 jet, and Ybarra and 
Lada| ( 2009 ) for HH 54. These data sets were also 
used in Paper I to investigate the nature of shocked 
emission observed in IRAC bands. The median of 
the total integration time per pixel ranges from 
21-322 s depending on the target. 

3. Image Analysis 
3.1. Method 

Based on Paper I, we expect the /4.5pm /-/s.o/xm 
ratio allows us to discriminate regions at high 
temperature or densities from the others since 
the 4.5-/im emission is more enhanced than 
the 8.0-/xm emission. In other words, a high 
^4.5/itm/^8.o/im ratio would be observed in regions 
at high temperature or density. However, this flux 
ratio can suffer severely from diffuse extended 



PAH emission at 8.0 fj,m (e.g., Reach et al. 2006 
Looney et al.|2007 ). When flux ratio maps are cal 



culated, masking the low signal to noise regions is 
standard practice, due to the relatively large un- 
certainties yielding very large (positive and nega- 
tive) ratios. For our purposes, this means that we 
could not easily compare the regions of shocked 
emission to the background regions in order to 
determine the origin of the flux ratio. 

Therefore, we alternatively apply the following 
approach: subtracting a flux-scaled 8.0-/im im- 
age from a 4.5-^m image (/4.5pm - r X Ig.o/im, 
where /4.5/im and Is.o/um are intensity at 4.5 and 



8.0 /im, respectively, and r is the scaling factor). 
The scaling factor r corresponds to the critical 
^4.5/im/4.o/iin ratio for discriminating between 
regions with high and low Ii.5fj,m/ h.Ofim'- the 
regions with h. 5 fj,m/ h.afim > r show Lt.s/im - 
r ^8.o/j,m values larger than the background; those 
with /4.5/xm/ h.afim < r show /4.5/im - rh.o^m 
values smaller than the background. The advan- 
tage of this method over the Zj.s^m/^s.o/xm map 
is that it is applicable even for regions with non- 
uniform diffuse background emission, in particular 
for HH 212 among our objects; and some regions 
in the L 1157 jet. This is because (1) such maps 
allow us to display the background regions with- 
out S/N masks; and (2) the measurements can be 
made relative to the value for the adjacent back- 
ground. This is in principle the same approach 
as that used for decades to investigate ionization 
and inferred shock conditions in Herbig-Haro ob- 



jects observed in atomic and ionic lines (Rcipurth 


and Heathcote||1992 Heathcote et al. 1996 


1998 


Fridlund and Liseau| 1998| |Hartigan et al. 


2000 


[Reipurth et al.||2002 


McCoey et al.||2004 


)■ 



Figure [T] shows an example of analysis in the 
HH 54 bow-shock region, which does not clearly 
show the presence of diffuse PAH emission in any 
band. Before obtaining the lis^m/^.o^m and 
/4.5pm — rls.o^m maps, the median value of the 
background is subtracted from each image, and 
the 4.5- and 8.0-/im images are convolved with 
the PRF at 8.0 and 4.5 /im, respectively. The 
PRFs at these two wavelengths are significantly 
different due to different diffraction patterns and 
degree of scattered light in the detector, and this 
approach is in particular useful for accurately com- 
paring their fluxes at the same angular resolution 
(FWHM—3" after convolution). Figure [I] shows 
that the /4.5/zm — ^-/s.o/im images allows us to suc- 
cessfully identify regions of high h.5[j,m/ h.Ofim ra- 
tios. Applying different r's (i.e., different critical 
h.5fj,m/ h.o/im ratios to show positive and nega- 
tive values) merely changes the contrast for dis- 
playing the features in shocks, as shown in later 
sections for the other regions as well. 

In the I4.5fj.rn — r ^8.o/zm images the foreground 
and background stars tend to show positive values 
as do shocked regions at high temperatures or den- 
sities. We discriminate between these two kinds 
of features using (1) conventional three color im- 
ages (blue, green and red for 3.6, 4.5 and 8.0 fj,m, 
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respectively), and (2) hs/im/ h.5/j,m flux ratio 
maps. The stars tend to appear blue in the three- 
color image, and exhibit a large /3.6/Ltm/^4.5/nm 
flux ratio (Paper I). The example for the HH 54 
bow shock region is shown in Figure [l] In Ap- 
pendix A we identify all the foreground and back- 
ground stars in our regions of interest. 

3.2. Results 

Figure [5] shows our whole view of the L 1157 
jet, HH 212, and HH 54. Figure [3] shows close- 
up views and the h.^^m — r/g.o/im images of the 
sub-regions selected from L 1157 and HH 212 in 
Figure [2] As for HH 54 in Figure [TJ the me- 
dian value of the background is subtracted from 
each image, and the 4.5- and 8.0-/im images are 
convolved with the PRF at 8.0 and 4.5 /zm, re- 
spectively, before obtaining the /4.5^m — ^s.o/xm 
maps. We mark prominent features in shocks in 
which /4.5/tm — ^s.o^m is higher than the adja- 
cent background using a letter and a number. For 
each region, the scaling factor r (i.e., the critical 
^4.5/xm/^8.o/Ltm value) is arbitrarily adjusted for 
contrast between the emission regions with high 
densities/temperatures and the remaining regions. 
The different r's required in different regions re- 
sults from varying shock conditions, yielding dif- 
ferent average or median /4.5^m/4.o/im ratios. 
This is independent of the diffuse background, as 
quantitatively proven in §4. 

The i4.5^,m — f^.o/xm images, in particular for 
L 1157 A, show a number of faint and unmarked 
point-like features with I^^m — r/g.o/im higher 
than the adjacent background. These are likely 
due to faint stars whose signal-to-noise for l4.5/j,m 
is less than 10 (see Appendix). In the text we 
limit our discussion to features clearly identified 
as shocks with this criteria for signal-to-noise. 
Some regions in Figure [3] show global variation of 
^4.5/Ltm — rlg.o^m in the background region due to 
the diffuse PAH emission apparent in the is.o/xm 
band. 

The h.5/j,m-rIa.oiJ,m images in Figures [TJ and |3] 
exhibit a number of knotty features in shocks with 
I4.5iu.rn — f-^.o/im higher than the background. As 
discussed in §2, these are the regions with rel- 
atively high temperature or density. The peak 
positions of the knots are listed in Table [TJ As 
described below in detail, these are categorized 
into: (1) those bracketed by bow shocks or bubble- 



like structures with small h.^^m. — r/s.o/im values 
(i.e., emission at relatively low temperatures or 
densities) ; (2) those located at the head of bow 
shocks; and (3) the others. Most are associated 
with bright knotty structures in the three-color 
image. 

In the HH 54 A, L 1157 A and C regions, HH 
54 A1-A4 and L 1157 A1/A2/C1 are bracketed by 
bow shocks or bubble-like structures with small 



'4.5/im 



rh.oytm values, i.e. 



emission at rela- 



tively low temperatures or densities. In HH 54 
A, four prominent knots are identified in two bow 
shocks which overlap. HH 54 Al and A3 are lo- 
cated at the head of two bow shocks. The po- 
sitions of HH 54 A1-A4 match peaks in the 4.5- 



/um. image identified by Ybarra and Lada (2009) 
(I, B, E, A, respectively, in their Figure 3.) In 
the L 1157 A region, multiple knotty structures 
are observed in Al and A3, and each of Al and 
A3 is bracketed by an asymmetric bow-shock-like 
structure at low temperature or density. In the 
L 1157 C region a single knot with a high tem- 
perature or density (CI) is bracketed by a single 
bubble-like emission feature at low temperature or 
density. Among the warm/dense knots described 
above, HH 54 A2, A4, L 1157 A3 and CI exhibit 
corresponding bright yellow knotty structures in 
the three-color images. 

The three-color image of the the HH 212 A 
and B regions exhibits bow shocks with angular 
scales of ^30", and the warm/dense knots Al and 
Bl are respectively associated with their heads. 
Another two warm/dense knots (A2, B2) are ob- 
served in the west wake of the HH 212 A/B bow 
shocks, respectively. The remaining part of the 
wakes, including the east wake of the HH 212 A/B 
bow shocks, show h^/jm — rIs.oiu.rn values smaller 
than the background, indicating that these regions 
have relatively low ^4.5/im/^s.o^m and thereby low 
temperatures or densities. The morphology in the 
upper stream of A2 and B2 suggest they are the 
heads of asymmetric small bow shocks (see also 
the H 2 1-0 S(l) high resolution images obtained 
by, e.g., |Zinnecker et al.|1998||Takami et al.|2006 |. 

Although Figure [3] is useful for identifying the 
remaining warm/dense knots (L 1157 A3/A4/B1- 
B5/C2 and HH 212 C1-C3), it does not clearly 
show how these are associated with gas of lower 
temperatures or densities. Figure [4] shows a close- 
up view of L 1157 A3, Bl, B3, B5, C2, and HH 
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212 C3. This figure shows that these warm/dense 
knots are associated with tails of smaller i4.5^m — 
r/g ounij thereby of lower temperature or density, 
towards the upper stream. 

Figure [5] shows the distribution of warm/dense 
knots over the entire region of the L 1157 jet and 
HH 212. The distribution of warm/dense knots 
is asymmetric between the northern and southern 
lobes of the L 1157 jet, as expected from previous 
observations of this object (e.g., Bachiller et al. 
2001 Looney et al.||200"7 ). A number of knots are 
observed in the northern lobe, distributed within 
an opening angle of about 30°. The regions with 
relatively low temperatures or densities, which 
look bright in the Ii.z^m — rls.ofim image, show a 
relatively continuous morphology analogous to the 
CO J =2 1 emission observed by Bachiller et al 
(2001). In contrast, only two bright warm/dense 



knots are identified in the southern lobe. The 
bubble-like morphology is not clearly observed in 



the CO J=2-l image by |Bachiller et al.| (|2001|) 



but is apparent in the CO J=l-0 image by|Gueth 



et al. ( 1996 ) at an angular resolution similar to 



our study. 

In contrast to the L 1157 jet, the warm/dense 
knots in HH 212 show a simpler distribution in 
Figure[5] Al, Bl, CI, C3 are distributed within 1° 
of a specific P.A. (24/204°), indicating a straight 
jet axis. In addition to the known positional sym- 



metry between Al and Bl (NK3 and SK3 in Zin- 



necker et al. 1998), A2 and B2 show an intrigu- 
ing symmetry with respect to the jet axis. Both 
are located at similar distances from the protostar 
(72" and 65", respectively) and with similar off- 
set angles from the jet axis (7°), but A2 is offset 
clockwise while B2 counterclockwise. 

4. Spectral Energy Distributions and Com- 
parisons with IRS Spectra 

To investigate the different physical conditions 
between warm/dense knots and cool/diffuse re- 
gions, we measured the fluxes in the four IRAC 
bands at the positions shown in Figures [6}|8] 
These include all the warm/dense knots we identi- 
fied in Figures [T] and [3] except HH 212-C1, in which 
the shocked emission is severely contaminated by 
non-uniform background emission at 5.8 and 8.0 
fim. For each position, we carefully subtracted 
the diffuse background emission by measuring it 



in the x- and y-directions and fitting it using a lin- 
ear function. We then measured the average flux 
density in a circular 6" aperture. Table [2] shows 
the measured fluxes at these positions. The un- 
certainty shown in Tableware due to non-uniform 
background emission in addition to the shot noise 
and readout noise of the detector. According to 
the IRAC Data Handbook 3.0, the absolute flux 
of extended emission measured using this camera 
is highly uncertain, and this depends on the spa- 
tial distribution of the emission both in and out- 
side the aperture. This fact is not included in the 
uncertainties listed in Table [2] One would expect 
those additional uncertainties to be ~ 20 % for 
the 3.6- and 4.5-/xm fluxes; ~ 30 % for 5.8-/xm; 
- 40 % for 8.0-/xm; - 0.3 mag. for the [4.5]-[8.0] 
color; and ~ 0.2 mag. for [5.8]-[8.0], based on the 
flux correction factor for different apertures listed 
in the IRAC Data Handbook. 

Figure [9] shows the spectral energy distributions 
(SEDs) at each position marked in the individual 
regions (L 1157 A/B/C, HH 54, HH 212 A/B/C). 
The SEDs for warm/dense positions show shal- 
lower slopes than those in cool/diffuse positions in 
most of the regions. Some or all of the SEDs for 
the warm/dense positions in L 1157 B, L 1157 C, 
HH 54, HH 212 A and HH 212 B show excess emis- 
sion at 4.5 /mi, while the other SEDs are relatively 
straight in Figure [9j In HH 212 C the differences 
in SED between warm/dense and cool/diffuse re- 
gions are marginal. 

These SEDs are similar to those measured in 
the six jet systems (including the L 1157 jet, HH 
54, and HH 212) in Paper I. In Paper I we showed 
that (1) the SEDs without 4.5-/mi excess can be 
explained well by shocked H2 emission with a 
power law cooling function (A oc T~~ a ); (2) the 
different SED slopes can be attributed to differ- 
ent power indexes or densities; and (3) those with 
4.5-/im excess require a contaminant at this wave- 
length, presumably CO, indicative of higher densi- 
ties than those without 4.5-/Ltm excess. To investi- 
gate this in detail, we plot the observed and mod- 
eled colors in color-color diagrams in Figures |10| 
and [TTj As in Paper I, models for H 2 emission are 
made for the following cases: (1) iso-thermal cases, 
and shock slabs with a power-law cooling function 
(A oc T~ a ); (2) thermal collisions with H+He and 
H2+He, corresponding to the cases with relatively 
high (^ 0.002 — 0.02) and low dissociation rates 
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(< 0.002-0.02, see Paper I). As modeled in Paper 
I, we adopt A-coeffcients provided by |Wolniewicz| 
et al. ( 1998 ), and collisional rate coefficients for H 2 



and He by Le Bourlot et al. ( 1999 ). For collisional 



rate coefficients with H, we adopt 



(2007) and Le Bourlot et al. (19991, and show 



Wrathmall et al. 



the results separately. According to [Wrathmall 
et al.l (|2007h , they provide rate coefficients with 



a better accuracy than Le Bourlot et al. (1999) 



due to improved representation of the vibration 
eigenf unctions. In contrast, coefficients provided 



by Le Bourlot et al. (19991 can explain the ob- 



served IRAC colors better in Paper I. The fluxes 
for H 2 +He collisions are calculated with the lowest 
36 energy levels (32 lines in the four IRAC bands), 
while those for H+He collisions are done with the 
lowest 49 levels (45 lines in the four IRAC bands) . 
The number of transitions we include is limited 
by the availability of collisional rate coefficients. 
The level populations of ortho- and para-H2 are 
calculated separately, and the resultant fluxes are 
combined assuming an ortho/para ratio of 3. The 
applicability and limitations of those calculations 
are discussed in Paper I in detail. 

Figures [T0| and [TT| show that, as in Paper I, the 
models with a power-law cooling function explain 
the observed colors at the majority of the above 
positions. In Figure [ITJ they are located near the 
colors for local thermal equilibrium (LTE) with a 
power index of the cooling function a of 3-6. The 
colors measured at most of the remaining positions 
show [4.8]-[8.0] smaller than the models, which is 
presumably due to contaminating emission from 
CO in the 4.5-ptm band (Paper I). As shown by 
Neufeld and Yuan| p008| , the contribution of CO 
to the entire 4.5-/zm flux is significant at densities 
> 10 7 cm~ 3 . The [4.5]-[8.0] colors larger than 
the LTE positions, indicative of lower densities, 
are observed at a few cool/diffuse positions in HH 
54. Figure [10] shows that the colors measured at 
all the positions have [4.5]-[8.0] smaller than the 
isothermal cases by at least 1-2 mag. 

In Figure [HI the [5.8]-[8.0] and [4.5]-[8.0] col- 



ors measured in each position show a correlation 
similar to the models with different power law in- 
dexes a. In particular, the warm/dense positions 
show smaller [5. 8] -[8.0] and [4. 5] -[8.0] colors, cor- 
responding to small a. This can be explained 
if the warm/dense knots include a larger frac- 
tion of warm gas that cool / diffuse regions shown 



in the h.sfim/ Is.Ofim images. In the HH 54 re- 
gion, the colors at the warm/dense positions in- 
dicate that the gas density is sufficiently high for 
LTE, while those at the cool/diffuse positions in- 
dicate lower densities 



we 



In Figures [10] and 1 1 
also show the /4.5ftm/4.o/im ratios correspond- 
ing to the scaling factor (critical flux ratio) r in 
Figures [Jj and [3j In each region, the [4.5]-[8.0] 
color corresponding to the middle r is larger than 
those at warm/dense positions, and smaller than 
those at cool/diffuse positions. This implies that 
the /4.5/xm/^8.o/zm ratios at the warm/dense po- 
sitions are larger than the middle r, while those 
at cool/diffuse regions are lower. All the results 
described above are consistent with our argument 
in §3.1 that /4.5/xm — ?^8.o/j,m allows us to discrim- 
inate between these warm/dense and cool/diffuse 
regions, independent of diffuse background at 8.0 
/Ltm. 



Neufeld et al. 


( 


2009 


) and 


Nisini et al. 


( 


2010) 



have made Spitzer IRS observations of the L 1157 
jet at the 5.2-37 /im range, and analyze the phys- 
ical conditions using H2 S(0) to S(7) lines. In Ta- 
ble [3] we compare our results in the same outflow 
with theirs in the following manner: (1) compar- 
ing the LTE temperature inferred by the H2 0-0 
S(5)/S(6)/S(7) lines ( jNisini et al.[[2010| with that 



obtained from the observed /; 



5.8/xm/^8.o/im 



ratio; 



(2) estimating the IRAC 8.0-/zm flux assuming 
that the flux is dominated by H 2 0-0 S(4) and 
S(5) lines ([Neufeld and Yuan|[2008j |Neufeld et al 



2009 ) and comparing with the observed IRAC flux 
at the same position; and (3) estimating the IRAC 
4.5-^im flux using the 0-0 S(9)/S(5) flux ratio, as- 
suming that the 0-0 S(9) is responsible for all the 
IRAC 4.5-/xm flux, and in LTE with the tempera- 
ture provided by Nisini et al. ( 2010[ ). The positions 
for the comparisons listed in Table [3] are selected 
based on signal-to-noise for the IRS spectra for the 
measurement of temperature in Figure 5 of Nisini 

eFaTT| poTol ) . 



In Table [3] the temperatures obtained from the 
IRS spectra and IRAC colors agree at most of 
the positions. In each region both temperatures 
are higher for the warn/dense knots than those 
for cool/diffuse regions we identified in Section 3. 
Furthermore, the IRAC 8.0-fim fluxes estimated 
from the 0-0 S(4) and S(5) lines approximately 
agree with the observed fluxes listed in Table [3] 
These support our discussion regarding modeled 
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colors which assumes that the emission at 5.8 
and 8.0- /im is dominated by H2, in particularly 
the 0-0 S(4) to S(7) lines (cf., |Smith and Rosen 



and this can be explained if the ejecta from the 
gas has a higher density than the ambient gas, 
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Ybarra and Lada 2009). In contrast, the IRAC 
4.5-^tm fluxes estimated in the above manner are 
systematically lower than the observed fluxes, by 
a factor of 0.06-0.3. This discrepancy can be at- 
tributed to a combination of the following for the 
estimated flux: (1) assumption of iso-temperature, 
which yields systematically large [4.5]-[8.0] col- 
ors (thereby low Ii.5/j,m/ la.o/jm flux ratios) than 
shock models with a power-law cooling function 
(Figures 10 and 11); (2) the fact that the other 
H2 lines and CO emission are not included. 

5. Discussion 

As described in §3 in detail, most of the 
warm/dense knots shown in Figures [1} [3] and 
[4]are categorized into: (1) those bracketed by bow 
shocks or bubble-like structures seen in regions 
at relatively low temperatures or densities ; (2) 
those located at the head of bow shocks; and (3) 
those associated with cometary tails. Throughout 
§5.1, we suggest that these are located at work- 
ing surfaces where the ejecta from the protostar 
interact with the ambient gas. In §5.2 we discuss 
the implications for the distribution of the ejecta 
from the protostar. 

5.1. Interaction between the jet and am- 
bient gas 

As shown in §3, several of the warm/dense 
knots we identified (HH 54 A1-A4, L 1157 Al, 
A2, CI) are bracketed by a bow shock or bubble- 
like feature at lower temperature or density. A 
likely explanation for these warm/dense knots is 
that these are reverse shocks into the ejecta from 
the protostar (i.e., Mach disks or cloudlet shocks 
Hartigan|1989 for schematic views) 



see, e.g. 



Such shocks have been observed towards a few pro- 
tostellar jets seen in atomic or ionic emission (see 



Reipurth and Bally||2001 for a review). 




Reipurth and Heathcote 


(1992 


); Fridlund and 


Liseau 


(1998 


); Reipurth et al. 


(2002|) have shown 



different excitations between Mach disks and for- 
ward bow shocks using the Ha and [S11] lines. In 
HH 34, Ha is more enhanced in the bow shock 
while [S11] is more enhanced in the Mach disks, 



yield 


ing lower ionization ( Reipurth and Heathcote 


1992 




Reipurth et al. 


2002 


). In contrast, the op- 



posite trend is observed in the L 1551 bow shock 
( Fridlund and Liseaup 998). The higher tempera- 
tures or densities in HH 54 A1-A4, L 1157 Al, A2, 
and CI coupled with the bow shock or bubble- like 
feature can be attributed to the different pre-shock 
densities of the ejecta from the protostar and the 
ambient gas. 

In HH 212 A and B the emission at high tem- 
peratures or densities is associated with the heads 
of bow shocks, while that at low temperature or 
density is associated with its wake (Figure[3| . This 
agrees with previous studies of bow shocks. Since 
we expect higher shock velocities at the head, this 
gas should have higher temperatures (and perhaps 
higher densities) at these positions ( Hartigan et al. 
1987||Reipurth and Bally|2001| for a review). The 
structures seen in HH 212 A and B suggest the 
presence of a relatively small bow shock (A2, B2) 
in the wake of a larger bow shock (Al, Bl). The 
structure of the small bow shocks is not clear at 
the angular resolution of our study, but is seen 
more clearly in the H2 2.12 /im images at higher 
angular resolutions ( Zinnecker et al.|1998 Takami 
eTaLl[2006l ). 

None of the bow shocks associated with HH 
212 A1/A2/B1/B2 clearly show evidence for the 
presence of reverse shocks, in contrast to HH 54 
A1-A4, L 1157 Al, A2, and CI. This may be at- 
tributed to a number of reasons. First, the reverse 
shocks can be significantly fainter if their preshock 
density is significantly lower than that of the am- 
bient gas (by a factor of more than 10-100 — see 
Hartigan 1989J. Second, reverse shocks may be 
relatively unstable, and thereby highly variable, as 



observed in HH 47 in optical emission lines (Har- 
tigan et al.||2005[ ). 



Cometary structures like those associated with 
L 1157 A3/B1/B3/B5, L 1157 C2, and HH 212 C3 
are rarely seen in protostellar jets. Indeed, none 
of the reviews of protostellar jets over the past 



decade has reported them (Hartigan et al. 


2000 


Eisloffel et al. 


2000 


Reipurth and Bally 


2001 


Bally et al. 


2007 


Ray et al. 


2007 


). Even so, it 



is natural to interpret these as ejecta interacting 
with the ambient gas at these positions, provid- 
ing a higher temperature or density and thereby 
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enhancing the 4.5-/im emission. The tails, which 
have a temperature or density lower than the head, 
could be gas entrained due to this interaction. 

5.2. Implications for jet driving 

According to the present jet paradigm, the 
ejecta of the jet is distributed as either a con- 
tinuous flow or a chain of knotty structures, mov- 



ing in opposite directions (e.g. Shu et al.| 


2000 


Konigl and Pudritz 


2000| |Hayashi et al.| 


1996 



Goodson et al.||1999[ ). As already stated, the 
observations of protostellar jets show a much 
more complex variety of features (knots, bow 
shocks, bubbles etc.) than the idealized mod- 
els. These features are often attributed to hy- 
drodynamic / magnetohydrodynamic instabilities 
of the shocked flows (e.g., Vishniac 1983 Wardle 
1990 Stone et aL||1995 ) , or time variability of the 
ejection in mass and direction (e.g., Raga et al 
2004b|a[[2W7al[2009| . 



Some warm/dense knots show distributions 
similar to the fragmented structures seen in hy- 
drodynamical simulations. The Kelvin-Helmholz 
and/or Rayleigh-Talor instabilities can fragment 
the ejecta along and across the flow axis (e.g., 



Ouyed et al. 2003 Klein et al. 1994), and these 



could explain fragmented structures in each of HH 
54 A, L 1157 A1-A3 and HH 212 C3. Perhaps 
cometary structures in the L 1157 jet could be at- 
tributed to fragmentation similar to that seen in 
wind (e.g., Cunningham et al. 2005), bullet (e.g., 



Raga et al. 



2007b), or collimated pulsed jet with 



precession (e.g., Rosen and Smith 2004 Smith and 



Rosen 20071. In contrast, none of these simula- 
tions show the following observed features. These 
are:- (1) two bow shocks nearly overlapping with 
each other in HH 54; (2) bubble-like and cometary 
features parallel to each other in the L 1157 C re- 
gion; (3) symmetry between HH 212 A2 and B2 
about the jet axis. Furthermore, it is not clear 
why the same possible instability in the L 1157 
B region does not break up the flow structures 
further downstream (i.e., the L 1157 A region). 

The symmetry between HH 212 A2 and B2 
about the jet axis is particularly notable, since the 
high degree of spatial symmetry along the jet axis 
of HH 212 also appear to rule out the possibil- 
ity that the features are imposed on an originally 



uniform jet via flow instabilities (Zinnecker et al. 



1998 ) . Zinnecker et al. ( 1998 1 pointed out that it is 



much more likely that those features arise through 
time variability at the driving source. Similarly, 
the symmetry between HH 212 A2 and B2 can be 
easily explained if it is due to the driving source. 

The above features in HH 54 A, L 1157 C and 
HH 212 A/B would therefore require an alterna- 
tive explanation. One possibility is that the mass 
ejection from these protostars is better approxi- 
mated by ballistic bullets ejected in roughly (but 
not exactly) opposite directions. Indeed, such a 
scenario can easily explain the observed structures 
described above. Although the detailed mecha- 
nism for such mass ejection is not clear, it has 
been proposed for the high mass protostellar out- 



flow OMC-1 ( e.g., |Allen and Burton||1993| |Nis- 
sen et al.|2007 l and proto-planetary nebulae (e.g., 
Dennis et al.|2008 ) . In the case of low-mass proto- 
stellar jets, Hartigan et al.| ( |2005[ ) revealed a num- 
ber of knotty structures in the well-studied bow 
shock HH 47 A, and argue that these are explained 
as small bullets that pass from the jet through 
the Mach disk and working surface to emerge as 



bumps in the bow shock. Yirak et al. ( 2009 ) have 



recently questioned the stability of a continuous 
magnetized fluid, and alternatively simulates the 
jet using a model similar to the one we propose 
above. 

If the trajectories of the bullets are approxi- 
mately the same as the arrows shown in Figure |4j 
the faintness of some features (L 1157 A3 and B5) 
can be explained as follows. The ambient gas is 
swept up by shocks associated with Al and A2. 
The passage of these shocks then allows the ambi- 
ent gas to move downstream, thereby weakening 
the impact of A3 and B5. Throughout, these fea- 
tures are consistent with the scenario of the shot- 
gun-like mass ejection discussed above. 

6. Conclusions 

We have analyzed data for three Herbig-Haro 
complexes (the L 1157 jet, HH 212, HH 54) ob- 
tained using the Infrared Array Camera (IRAC) 
on the Spitzer Space Telescope to investigate how 
these flows interact with the ambient gas. A 
combination of 3.6-, 4.5- and 8.0-fim images has 
allowed us to successfully identify a number of 
warm/dense knots in shocked emission in molecu- 
lar gas at 1000-4000 K. Most of them are catego- 
rized into: (1) those bracketed by bow shocks or 
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bubble-like structures seen in regions at relatively 
low temperatures or densities; (2) those located 
at the head of bow shocks; and (3) those associ- 
ated with cometary tails. Our results show that 
these locations are distributed not only along the 
jet axis, as expected, but also across it. 

It is likely that the first group of warm/dense 
knots are associated with the ejecta from the 
protostar (i.e., Mach disks or cloudlet shocks). 
Analogous to some optically-visible Herbig-Haro 
objects, different h.sfxm/Is.Ofim ratios between 
the reverse and forward shocks can be attributed 
to the different conditions of the pre-shock gas, 
i.e., the ejecta from the protostar and ambient 
gas. The third group is explained if the ejecta 
interacts with the ambient gas at the location 
of the warm/dense knots, providing a relatively 
high shock velocity and thus compressing the gas. 
Throughout, we suggest that all of the above fea- 
tures are located at the working surfaces where the 
ejecta from the protostar interact with the ambi- 
ent gas. 

The origin of some knotty structures can be 
attributed to shock instabilities previously mod- 
eled by other authors. In contrast, others do not 
seem to be easily attributed to the same physi- 
cal processes with a continuous flow or a chain of 
well aligned knotty structures. A simple alterna- 
tive explanation for these observations, in partic- 
ular for the HH 54 A, L 1157 C and HH 212 A- 
B regions, is that the mass ejection is more akin 
to a "shot gun" , periodically ejecting bullets of 
material along similar but not identical trajecto- 
ries. This explanation challenges to some degree 
the present paradigm for protostellar jet flows as 
quasi-continuous, but give clues to reconciling our 
understanding of the mass ejection mechanism in 
high and low mass protostars and evolved stars. 

We are grateful to Dr. D. Ncufcld for use 
of their IRS spectra, and to Dr. F. Shu, R. 
Krasnopolsky, C.-F. Lee, O. Morata, M. Cemeljic, 
N. Hirano, and S. Takakuwa for useful discus- 
sions. We also thank the anonymous referee 
for useful comments. The IRAC images were 
obtained through the Spitzer Archive operated 
by IPAC. This research made use of the SIM- 
BAD data base operated at CDS, Strasbourg, 
France, and dNASA's Astrophysics Data System 
Abstract Service. This study is supported from 



National Science Council of Taiwan (Grant No. 
97WIA0 100327). 

Facilities: Spitzer Space Telescope (IRAC). 
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A. Identification of foreground and background stars 



Figure [T2| shows the foreground and background stars identified in the L 1157 A, B, and C regions. As 
described in the caption of Figure [T] identification of these stars is based on their point-source appearance 
and at least one of the following: (1) a blue color in the three-color image, and (2) a higher /3.6/xm/^4.5/zm 
flux ratio than the surrounding region. The last is based on the detailed study of shocked emission in 
Herbig-Haro objects of Paper I 0$0.6 for Herbig-Haro knots; ~1.6 for foreground stars with Ak—0; ^0.9 
for background stars with and here we apply this to the regions where the signal-to-noise is larger 

than 10 for I^.^^m- Contaminating emission from shocks can yield a lower /3.6/Ltm/^4.5/zm flux ratio for some 
stars, hence we apply a more qualitative identification with the other criteria. 

For HH 212, the 3.6-/im image suffers significantly from non-uniform diffuse extended emission, which 
would affect identification of foreground and background stars using the Is.e/jun/Ii.s/im flux ratio map. This 
diffuse emission component is also seen at 8.0 /jm with a very similar distribution, and can be removed 
reasonably well if we scale the flux of the 8.0-/im image and subtract it from the 3.6-^tm image. 

Figure 13 shows how this process works: i.e., the image at 3.6 /im, that at 8.0 /Ltm the flux of which is 
scaled to a level similar to 3.6 /im, and the 3.6-/mi image from which the scaled 8.0-/im image is subtracted. 
Foreground and background stars are faint in the scaled 8.0-fim image as expected: they should be 2% and 4% 
of the 3.6-^m flux for most spectral types (i.e., [3.6]-[8.0]~0) and extinction A K =Q and 5, respectively, based 
on Chapman et al. (2009). In contrast, shocked emission is apparent at 8.0 /im. This implies our imaging 
process keeps the stellar fluxes high, while it degrades shocked emission. Throughout, this would allow an 
even clearer discrimination between shocks and foreground/background stars using the hs/im/ h.5fj,m flux 
ratio (the exact formula for this case is (/3.6/im — ax /8.o/im)/^4.5/im, where a is the scaling factor for the 
8.0-^m flux). Figure 14 shows the foreground and background stars identified in the HH 212 A, B, and C 
regions, replacing the h.6fj,m/ h.5fim images with (I 3 . 6 fim ~ a x h.Ofim) / h.5^im images. The scaling factor a 
was determined using the least squares method to minimize the residual of subtraction for the non-uniform 
diffuse extended emission. 

The most likely origin of the non-uniform diffuse extended emission is PAHs illuminated at the surface 
of the molecular cloud. The measured scaling factor a of 0.071 (HH 212 A) and 0.0821 (HH 212 B,C), 
which correspond to the /3. 6/xm/ h.oum, is in excellent agreement with recent models for UV-excited PAH 
emission by |Draine aiiaTrj pfJTj7] (0.063-0.084 for {7=1-10 and 0.059-0.083 for [7=100, where U is the UV 



field normalized to the local interstellar radiation field). 

All the stars identified in Figure [l] are seen in the 2MASS images for HH 54, for which the images with a 
long exposure (6x) are available. Most of the stars in Figure 12 and some of the stars in Figure 14 are seen 



in the 2MASS images, while the remaining faint stars are not, due to the significantly shallower detection 
limit of 2MASS. 



1 This column density is higher than the maximum column density measured in some molecular clouds using background stars 
l |Chapman et al.|2009[ ). See also Mot te and Andrij ( |200l| for the molecular envelopes of low-mass protostars. 
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Fig. 1. — The HH 54 bow shock region, (top-left) three-color image (blue, green, red for 3.6, 4.5 and 8.0 
/um, respectively) (top-middle) the 4.6/im/^4.5/im flux ratio map. The ratio is shown in regions where the 
signal-to-noise is larger than 10 for Z^.s^rn- (top-right) A.s/im/^s.o/Ltm (bottom) l4.5fj,m. — rls.0fj,m images with 
different r's. The warm or dense knots in shocks (i.e. those with large /4.5/xm — ^s.o/zm and -?4.5^tm/-f8.o/xm 
values) prominent in the i4.5^tm — rIs.o/j,m and Ji.s^rn/ h.o/im images are marked with dashed circles and 
ellipses. Stars identified based on the /3.6/^m/^4.5/im ratio and the three-color image (see §2 for details) are 
marked with dotted circles in the /3.6/im/^4.5^m map, and white asterisks in the /4.5^m — r ^8.o/xm maps. 
Dashed curves in the three-color image indicate two bow shocks discussed in Sections 3 and 4. 
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Fig. 2. — Three-color images (blue, green, red for 3.6, 4.5 and 8.0 /xm, respectively) for the L 1157 jet, HH 
54 and HH 212. Dashed boxes show the regions where we apply the analysis described in §2. Green crosses 



in the L 1157 and HH 212 regions indicate the position of the protostar, based on Bachillcr et al. (2001) 



and Lee et al. (2008), respectively. 
Garatti et al.||2006[ ) 



The driving source of HH 54 has not been clearly identified ( Caratti o 
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15 



Fig. 3. — Three-color images (blue, green, red for 
3.6, 4.5 and 8.0 /im, respectively) and /4.5^rn — 
r/g o^im maps with different r for the regions se- 
lected in Figure[2](for L 1157, HH 212). All images 
are oriented to make left the downstream direc- 
tion. Dashed ellipses show the warm/dense knots 
in shocks prominent in the /4.5/xm — ^^8.o/im maps. 
The corresponding regions are also marked in the 
three-color image. These are labeled with a com- 
bination of a letter corresponding to the region 
(A-C), and a number (1-5). L 1157 A4 and A5 
are identical with Bl and B2, respectively, but are 
shown with different r values for the maps of L 
1157 A and B. The asterisks show stars identified 
using the three-color image and hs^m/h.s/im 
flux ratio (see §2 and Appendix). 
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20" 




Fig. 4. — Same as Figure [3] but for some 
warm/dense knots associated with tails (L 1157 
A3/B1/B3/B5/C2, HH 212 C3). For most of the 
regions the value of r has been changed to show the 
contrast between the head and tail(s) of cometary 
structures. 
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Fig. 5. — (a) Three-color image of L 1157. The arrows indicate the positiorTbf warm/dense knots marked 
in Figure [3] relative to the position of the protostar ( Bachiller et al.|[2001[ ) marked with the green cross, (b) 
The /4.5(im — rJg.o^m images in Figure [3| are superimposed on the 4.5-^m image in gray scale of the same 
region. The warm/dense knots are labeled in the same manner as Figure [3j and the arrows are the same as 
(a). The asterisks show the stars. (c)(d) Same as (a)(b) but for HH 212. The position of the protostar is 
based on Lee et al. ( 2008 ) . 
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20 



Fig. 6. — Positions for circular aperture photome- 
try (6") in the L 1157 jet, on the h.5nm — rls.0fj,m 
and 3.6-, 4.5-, 5.8-, and 8.0-^m images. The white 
and black circles are the apertures for warm/dense 
and cool/diffuse components, respectively. These 
are labeled as follows: WD — warm/dense posi- 
tion shown in the h^^m — r/s.o/im image; CD 
- cool/diffuse position shown in the /4.5^m — 
rh.Ofim image; A1/A2/B1-B5/C1/C2 — posi- 
tions (or adjacent cool/diffuse positions) for the 
warm/dense knots identified in Figure [3] Dashed 
boxes and lines show positional coincidence be- 
tween the i4.5^im — r/g.o/im and 3.6-^m images. 
Stars also show values larger than the surround- 
ing region in the /4.5/im — rIs.o/j,m, and these are 
marked in these images using asterisks (see §2 and 
Appendix A for identification). 
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Fig. 7. — Same as Figures [6] but for HH 54. A1-A4 in the labels are based on identification of warm/dense 
knots in Figure [TJ See Appendix A for identification of stars in the Ii.^^m — r/s.Ofim image. 
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/4.5nm — f"/8.»nm 




Fig. 8.— Same as Figure [6] but for HH 212. 
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Fig. 9. — Spectral energy distributions for the four IRAC bands at the positions shown in Figures [6]-[8] Solid 
and gray lines are those measured in the warm/dense knots (WD) and cool/diffuse regions (CD), respectively. 
Error bars are shown only for those larger than the size of the dots. These are based on Table [2] and possible 
systematic errors for the absolute flux calibration (see text) are not included. 
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Fig. 10. — Color-color diagram for observed fluxes 
in the individual regions, and isothermal H2 emis- 
sion. Solid and gray dots indicate the colors mea- 
sured in warm/dense and cool/diffuse regions, re- 
spectively. Error bars are shown only for those 
larger than the size of the dots. These are based on 
Table [2j and possible systematic errors for the ab- 
solute flux calibration (see text) are not included. 
Models are shown for the following three cases 
from left to right: (left) thermal collisions with H2 
and He, adopting the collisional rate coefficients 
provided by Le Boutlot et al. (1999); (middle) 
same but collisions with H and He; (right) same 
as the middle but adopting the rate coefficients 
provided by Wrathmall et al. (2007) for collisions 
with H. The [4.5]-[8.0] colors corresponding to r 
(i.e., the critical /4.5/xm/^s.o/mi ratios) in Figures 
[T] and [3] are shown at the right side of each row. 
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Fig. 11. — Same as Figure 10 but shock mod- 
els with a power-law cooling function (A oc T~ Q ) 
are plotted instead of using the isothermal as- 
sumption. The numbers shown in the plots 
(3.0/3.8/4.7/6.0) are the power index a. 
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. 12. — The three-color image, the -f3.6^m/^4.5/mi and the I^.s^m — ^s.o/mi maps for L 1157 A-C. The 
/im/^4.5/jm ratio is shown in regions where the signal-to-noise is larger than 10 for i^s^tm- 
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Fig. 13. — Subtraction of the non-uniform diffuse extended emission at 3.6 /jm for part of the HH 212 jet. 
(top-left) the 3.6-/im image before subtraction, (top-middle) flux-scaled image at 8.0 fim. (top-right) the 
3.6-/^m image subtracted from the flux-scaled image at 8.0 /im. The three figures arc displayed with the 
same contrast level as shown by the scale bars, (bottom) one-dimensional flux distribution extracted from 
the images before and after subtraction (gray and solid lines, respectively). The dot-dashed line indicates 
the zero flux level. 
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Fig. 14. — Same as Figure [12] but for HH 212 A-C. Non-uniform diffuse extended emission at 3.6 /im is 
subtracted using the flux-scaled image at 8.0 /im. 
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Table 1 

Peak Position of Warm/Dense Knots 



Object 


Feature 


R.A. (2000) 


Dec. (2000) 


Distance to the 
protostar (arcsec) a 


P.A. from the 
protostar (deg.) a 


T 11 ^7 
Li 110 ( 


A 1 b 


ZU.oo.04.0 


+00.1)4.0 / 


I/O 


ooy 






ZU.oO.Oo.y 


+00.U4.0Z 


1 (\R 
100 


o4U 






ZU.oO.OO.z 


+OO.U4.00 


1 fiQ 

loo 


o41 






90-38-57 a 


4-68-04-48 


1 


343 




A2 b 


20:39:00.1 


+68:04:21 


130 


345 






20:39:00.8 


+68:04:14 


122 


346 




A3 b 


20:38:58.9 


+68:04:13 


124 


341 






20:38:58.4 


+68:04:06 


118 


338 






20:38:58.4 


+68:03:58 


111 


337 




Bl (A4) 


20:38:57.6 


+68:03:32 


90 


328 




B2 (A5) 


20:39:00.5 


+68:03:34 


84 


338 




B3 


20:38:58.2 


+68:03:08 


68 


319 




B4 b 


20:39:00.6 


+68:03:17 


68 


333 






20:39:00.5 


+68:03:08 


61 


329 




B5 


20:39:04.9 


+68:02:55 


39 


350 




CI 


20:39:09.3 


+68:01:20 


59 


163 




C2 


20:39:12.1 


+68:01:11 


73 


153 


HH 212 


Al b 


05:43:53.9 


-01:01:29 


92 


24 






05:43:53.7 


-01:01:32 


88 


23 






05:43:53.8 


-01:01:35 


86 


25 




A2 


05-43-52 8 


-01:01:44 


72 


17 




Bl b 


05:43:49.0 


-01:04:09 


83 


205 






05:43:49.2 


-01:04:04 


78 


205 




B2 


05:43:49.2 


-01:03:49 


65 


211 




CI 


05:43:47.7 


-01:04:58 


137 


204 




C2 


05:43:48.9 


-01:04:36 


109 


200 




C3 b 


05:43:48.3 


-01:04:36 


113 


204 






05:43:48.4 


-01:04:33 


109 


204 






05:43:48.5 


-01:04:27 


103 


205 






05:43:48.7 


-01:04:26 


101 


203 


HH 54 


Al 


12:55:54.9 


-76:56:06 








A2 


12:55:51.0 


-76:56:21 








A3 


12:55:53.7 


-76:56:22 








A4 


12:55:49.5 


-76:56:30 







a We adopt (aj2ooo^j20oo)= (20:39:06.2,68:02:16) and (05:43:51.4 
driving sources, respectively ( |Bachiller et al. |2001| |Lee et aF 
known (Caratti o Garatti et al. 2006; Paper I) 



01:02:53) for the L 1157 jet and HH 212 
2008). The driving source of HH 54 is not 
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These features contains multiple peaks. 



Table 2 

Fluxes measured in warm/dense and cool/diffuse regions 



Region 


Position 




Flux (MJy str- 1 ) 








3.6 /jm 


4.5 /um 5.8 /im 


8.0 /xm 



L 1157 A 



L 1157 B 



L 1157 C 



HH 54 A 



HH 212 A 



HH 212 B 



HH 212 C 



WD-A1 

WD-A2 

CD-A1 

CD-A2a 

CD-A2b 

WD-B1 

WD-B2 

WD-B3 

WD-B4 

WD-B5 

CD-B2 

CD-B3B4 

WD-C1 
WD-C2 
CD 

WD-A1 

WD-A2 

WD-A3 

WD-A4 

CD-I 

CD-2 

CD-3 

WD-A1 

WD-A2 

CD-I 

CD-2 

CD-A2 

WD-Bla 

WD-Blb 

WD-B2 

CD-I 

CD-2 

CD-B2 

WD-C2 

WD-C3a 

WD-C3b 

CD-I 

CD-2 

CD-C3a 

CD-C3b 



0.115 
0.216 
0.054 
0.106 
0.047 

0.097 
0.386 
0.112 
0.847 
0.213 
0.077 
0.114 

2.383 
0.775 
0.668 

0.89 
1.21 
2.06 
1.53 
0.23 
0.37 
0.28 

0.25 
0.434 
0.06 
0.02 
0.17 

1.43 
1.48 
0.738 
0.08 
0.176 
0.108 

0.153 
0.12 
0.10 
0.085 
0.059 
0.079 
0.063 



0.005) 
0.005) 
0.005) 
0.005) 
0.005) 

0.011) 
0.011) 
0.008) 
0.008) 
0.004) 
0.011) 
0.008) 

0.009) 
0.009) 
0.009) 

0.04) 
0.04) 
0.04) 
0.04) 
0.04) 
0.04) 
0.04) 

0.04) 

0.010) 

0.02) 

0.02) 

0.02) 

0.03) 

0.03) 

0.008) 

0.02) 

0.009) 

0.008) 

0.012) 

0.05) 

0.05) 

0.006) 

0.015) 

0.013) 

0.013) 



0.006) 
0.005) 
0.006) 
0.005) 
0.005) 

0.007) 
0.007) 
0.010) 
0.010) 
0.007) 
0.007) 
0.010) 

0.013) 
0.013) 
0.013) 

0.03) 
0.03) 
0.03) 
0.03) 
0.03) 
0.03) 
0.03) 

0.03) 
0.02) 
0.013) 
0.013) 
0.04) 

0.06) 
0.06) 
0.02) 
0.04) 
0.015) 
0.013) 

0.306 (0.02) 
0.347 (0.04) 
0.302 (0.04) 
0.233 (0.011) 

0.13 (0.02) 
0.224 (0.016) 
0.184 (0.016) 



0.304 
0.488 
0.178 
0.298 
0.138 

0.186 
0.751 
0.207 
1.651 
0.474 
0.154 
0.258 

5.766 
1.768 
1.162 

2.26 
3.17 
4.05 
4.52 
0.63 
0.99 
0.64 

0.52 
0.87 
0.168 
0.081 
0.39 

4.00 
3.96 
2.02 
0.18 
0.388 
0.267 



0.59 
0.962 

0.56 
0.802 
0.422 

0.259 
1.163 
0.29 
2.21 
0.492 
0.353 
0.59 

4.85 
1.61 
1.72 

3.05 
4.94 
6.83 
3.84 
1.99 
2.79 
1.91 

0.32 
0.75 
0.24 
0.16 
0.46 

2.95 
2.85 
1.31 
0.30 
0.48 
0.42 

0.43 
0.49 
0.43 
0.44 
0.23 
0.36 
0.31 



0.03) 

0.009) 

0.03) 

0.009) 

0.009) 

0.012) 

0.012) 

0.02) 

0.02) 

0.010) 

0.012) 

0.02) 

0.03) 
0.03) 
0.03) 

0.09) 
0.09) 
0.09) 
0.09) 
0.09) 
0.09) 
0.09) 

0.03) 
0.02) 
0.03) 
0.03) 
0.05) 

0.07) 
0.07) 
0.04) 
0.04) 
0.04) 
0.02) 

0.03) 
0.05) 
0.05) 
0.02) 
0.03) 
0.04) 
0.04) 



0.95 
1.43 
1.20 
1.72 
0.93 

0.22 
1.22 
0.29 
2.17 
0.37 
0.73 
1.08 

4.53 
1.47 
2.49 

4.04 
4.39 
8.66 
3.17 
3.19 
4.33 
3.94 

0.25 
0.64 
0.40 
0.40 
0.52 

2.00 
1.70 
0.79 
0.31 
0.42 
0.43 

0.50 
0.66 
0.51 
0.62 
0.40 
0.64 
0.55 



0.03) 
0.02) 
0.03) 
0.02) 
0.02) 

0.02) 
0.02) 
0.03) 
0.03) 
0.02) 
0.02) 
0.03) 

0.05) 
0.05) 
0.05) 

0.04) 
0.04) 
0.04) 
0.04) 
0.04) 
0.04) 
0.04) 

0.05) 
0.04) 
0.04) 
0.04) 
0.08) 

0.07) 
0.07) 
0.03) 
0.05) 
0.07) 
0.02) 

0.03) 
0.07) 
0.07) 
0.02) 
0.07) 
0.07) 
0.07) 
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